UB-ECM-PF-03/27 
FIRENZE DFF-407/10/03 



Two tricritical lines from a Ginzburg-Landau expansion: 
application to the LOFF phase 

■ R. Casalbuoni la and G. Tonini 6 

O 

a Departament ECM, Facultat de Fisica, 
Universita de Barcelona and Institut de Fisica d'Altes Energies, 
O ■ Diagonal 647, E-08028 Barcelona, Spain 

b Dipartimento di Fisica, Universita di Firenze, 
O ', 1-50019 Firenze, Italia 

and 

I.N.F.N., Sezione di Firenze, 1-50019 Firenze, Italia 



> 
OO 



m 
O 



email: casalbuoni@fi.infn.it, tonini@fi.infn.it 
Abstract 



We study the behavior of the two plane waves configuration in the LOFF phase close to 
T = 0. The study is performed by using a Landau-Ginzburg expansion up to the eighth order in 
the gap. The general study of the corresponding grand potential shows, under the assumption 
that the eighth term in the expansion is strictly positive, the existence of two tricritical lines. 
This allows to understand the existence of a second tricritical point for two antipodal plane waves 
in the LOFF phase and justifies why the transition becomes second order at zero temperature. 
^ , The general analysis done in this paper can be applied to other cases. 



1 Introduction 



The condition for the BCS pairing is that the participating fermions should have opposite 
momenta. This condition is realized when the Fermi spheres of the two species are the same. 
However in various circumstances this does not happen. For instance, in ordinary matter, in the 
presence of a magnetic field the energies of the spin up and spin down electrons are different and 
this gives rise to a splitting of the two Fermi surfaces. In the case of quarks a mass difference 
between the two pairing species leads indeed to such a separation. The separation can also 
be produced by the condition of weak equilibrium and/or charge neutrality. All these different 
instances can be essentially described in terms of a difference in chemical potential (real or 
effective), 5u = — /Urf)/2, with \jl u and u,i the chemical potentials of the two fermions. As 
long as the separation is not too large the BCS pairing still arises. However increasing the 
separation leads to a first order phase transition to the normal phase at a point known as the 
Clogston-Chandrasekhar jTj limit, given by Sfi = Sfii = Aq/\/2, with Ao the BCS gap. 

1 On leave from Dipartimento di Fisica, Universita di Firenze. 
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Already many years ago Fulde, Ferrel [2] and Larkin, Ovchinnikov 3 pointed out that close 
to this first order transition a different pairing mechanism could arise 2 . In this case each of the 
fermions lie close to their own Fermi surface and therefore the pair has non zero momentum. 
This phase has been called both FFLO and LOFF in the literature. Here we will adopt the name 
LOFF 3 . As a consequence of the non zero momentum of the pair translational and rotational 
invariances are broken. It turns out that a possibility is the formation of crystalline structures 
with the vertices of the cell identified by the possible momenta of the pair Correspondingly 
the condensate has a spatial dependence of the type 

p 

A(r) = Y, A m e 2iq ™ r . (1) 

m=l 

The case of a single plane wave has been studied both in [5] and 0. In [3] two other cases 
were considered: the case of two antipodal plane waves and the one corresponding to a cube. 
All these studies have been performed at zero temperature, T = 0. More recently this problem 
has become of interest in QCD [7j just for the reasons mentioned above, that a splitting in 
the Fermi energies can be produced by different quark masses and/or by weak equilibrium and 
charge neutrality. In [S] 23 different crystalline structures have been studied (still at T = 0) with 
the conclusion that the most favorable structure would correspond to the face-centered cube, 
described by eight plane waves. 

As far as the studies at T ^ are concerned, the phase diagram in the plane (6/j,, T) 
illustrating the transition from the BCS phase to the normal phase can be found in j^j. The 
main result is the existence of a tricritical point, that is of a point where a first-order transition 
line and a second order one meet. This is easily understood since at 5/jl = we have the usual 
BCS second order transition at T c ~ 0.57Ao, whereas at T = we have a first order transition at 
Sfj, = Ao/y/2. The tricritical point, Ptric occurs at d^tric ~ O.6IA0 and T tr i c ~ 0.32Ao [TT^ [TT] . 
It turns out that this tricritical point is the end point also for any LOFF phase, in the sense 
that at that point also the vectors q m vanish and we are back to a BCS-like situation. A careful 
analysis of the LOFF phase has been done in The result is that it is very likely that the 
favored phase close to the tricritical point is the one corresponding to two antipodal vectors, 
that is to a condensate 

A(r) = 2Acos2q- r. (2) 

Furthermore the transition from this LOFF phase to the normal one is first-order close to the 
tricritical point. On the other hand, in [2] it has been shown that the transition is second-order 
at T = 0. 

From the previous considerations we see that there are at least two problems to understand, 
the first one is the possible change of the favored crystalline structure from the tricritical point 
to zero temperature and the second one is the change in the nature of the transition in the 
case of the antipodal structure. A partial answer to the first problem comes from two analysis 
made in the two-dimensional cases |121 IT3] . Essentially what these authors find is that there 
is a series of transitions starting from the antipodal case, close to the tricritical point, between 
different crystalline structures. In particular, in ref. ^3] it is argued that the complexity of the 
crystalline structure increases going toward T = 0. 

2 More recently other possibilities have been considered, see |1] 
3 For recent reviews of this phase, see |^JE] 
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In the present paper we will be concerned with the second problem. This has been studied 
in |14j using the method of quasi-classical Green's functions |15j . The result is that at very low 
temperature T ~ 0.043Aq ~ 0.075T C there is a further tricritical point allowing the transition 
from the first-order to the second order. In this paper we would like to study this problem by 
using the more classical approach based on the Ginzburg-Landau expansion around a second 
order phase transition. The complication here is that we have to expand the grand potential (or 
the free energy) up to the eight order in the constant A of eq. ©: 

Q = aA 2 + i/5A 4 + i 7 A 6 + ^A 8 . (3) 

In fact, whereas around Ptric the 7 coefficient is positive, and the expansion up to A 6 is enough to 
characterize the tricritical point, 7 becomes negative going down in temperature |14j . Therefore 
for reasons of stability of the theory a further positive term in the expansion is necessary. 
Unfortunately no calculation of the eighth order term exists so far, therefore we will perform our 
study assuming 5 positive up to zero temperature and we will express our results parametrically 
in 5 itself. 

Before studying the LOFF case we will consider the general phase structure for a grand 
potential f2 as given in eq. (jHJ). Since we assume 5 to be strictly positive, £1/5 depends only on 
the ratio of the three parameters a, (3 and 7 to 5. We will show that in this three-dimensional 
space there are two lines of tricritical points. The point Ptric hes on the line a = (3 = 0, 7 > 0. 
The other line is still in the plane a = but it lies on the 7 < side and it is given by the 
equation 9(35 - 2j 2 = 0. 

As we have already noticed, from ref. ^1] one can argue that a second order line should 
exist starting from T = and ending to another tricritical point. This point is located at low 
temperature T/Aq ~ 0.043. We will take advantage of this fact by evaluating the parameters 
a, f3 and 7 up to the second order in T/Ao- Then one has to map the space (a, /3, 7) on the 
parameter space of the theory, that is T, 5fi and the vector q characterizing the antipodal 
state to get the result. As already anticipated we find a second tricritical point with a location 
depending on the value of 5. We recover the results of ref. |14j for very small values of 5. In 
particular the temperature and the 5/i values characterizing this point vary of about 25% and 
1.5% in the range < 5 < 1. Therefore we confirm the results of the analysis made in |14j . as 
long as the parameter 5 is strictly positive. 

We would like to point out the analysis we have made of the critical points of the grand 
potential given in eq. © is completely general and as such it can be applied to all the situations 
where the eighth term in the Ginzburg-Landau expansion plays a role. 

In Section 2 we describe the Ginzburg-Landau expansion for the LOFF case of two antipodal 
plane waves. Section 3 is dedicated to the study of the critical points of the grand potential (|2J). 
In Section 4 we apply the previous considerations to the LOFF case. Conclusions are given in 
Section 5. 

2 The Ginzburg-Landau expansion for the case of two antipodal 

waves 

As discussed first in [Hj and then in |Hj (see also El ) one can derive the Ginzburg-Landau 
expansion for the grand potential by first expanding the gap equation in the gap parameter 
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and then integrating over the parameter itself. All the relevant equations are given in the 
literature at zero temperature. The extension at finite temperature is trivial and it is made 
by the introduction of the Matsubara frequencies. As already said we will be concerned only 
with the first three coefficients a, f3 and 7 which are the ones evaluated at T = in [3] and 
|B], and we will assume the coefficient 5 to be strictly positive. Notice that we are interested at 
the way in which the first order transition line coming from the tricritical point becomes second 
order close to zero temperature. Therefore it will be enough to evaluate the coefficients of the 
Landau-Ginzburg expansion close to T = 0. For that we will perform an expansion in powers 
of T/Aq where Ao is the BCS gap, that is the gap evaluated at a value of 5p smaller than the 
Clogston-Chandrasekhar value, dpi = Aq/\/2. 

The grand potential SI is given in the GL approximation up to the sixth order in the gap by 

1 P 

= "-( E I n (qk,q„) - 1] A£A n 5 qk _ qn 

fc,n=l 



+ \ E J (qk,q£,qm,q n )AfcA £ A^A n( 5 qk _ q , +qm _ qn 

k,£,m,n=l 
P 

+ 3 E ^(qk,q^qm,qj,qi,qn)A^A^A^A i A*A„5 qk _ q , +qm _ qj+qi _ qn ) , (4) 

k,£,m,j,i,n=l 

where P is the number of independent plane waves in the condensate and 

IK*.*) = + f / £ f_ S s « £°° § A MB, In. {<,}) , (5) 

«/(qi,q2,q 3 ,q4)= + J -£ J & d£ j |- JJ ftE, 6p, {q}) , (6) 



^(qi,q 2 ,q3,q4,q 5 ,q6)= + '-f J ^ f & J_ " ^ fl fi&, {q}), (7) 
where, by putting w = vpw, 



i=l 



ME, Sfi, {q}) = ■ : . (8) 

E + iesigaE _ §fi+{ _ m _ 2 ^ k= ^_ 1)kw . qk] 

Moreover the condition 

M 

£(-l)*«lk = (9) 

k=l 

holds, with M = 2,4,6 respectively for II, J and K. Furthermore p is the density of states at 
the Fermi surface given by p = gpp/ir 2 VF, with pp Fermi momentum, vp the Fermi velocity and 
g the number of degrees of freedom. In the case of two-flavor QCD in the color superconductive 
phase (2SC phase) one has g = 4, pp = p and vp = 1. If the vectors q define a crystalline 
structure they belong to the orbits of the point group of the crystal. In most of the cases 
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considered in [H] there is a single orbit, and this is the case for two antipodal waves. Therefore 
we assume that 

A k = A* k = A (for any k) . (10) 
Then we can rewrite (jlj as follows (for two antipodal waves): 

- = aA 2 + ^A 4 + J A 6 , (11) 
p 4 o 



where 



a = 2 i^iM, (12) 
99 



p 



2 , 

P = Yl J (qk,q£,qm,qn)<5q k -q,+q m -q n , (13) 



gp 



k,l,m,n=l 

2 



P 



-q^+qm-qj+qi-qn- 

(14) 

k,£,m,j,i,n=l 

These expressions, which are obtained in [S] for T = 0, can be easily extended to finite tempera- 
ture by transforming the integration over the energy into a sum over the Matsubara frequencies 
according to 

/+oo 
dE -> 2-KiT ^2 , E — > iuj n = (2n + l)vrT. (15) 
n=— oo 

The way of evaluating these integrals at T = is explained in [S]. One can follow the same 
method and perform the sum over the Matsubara frequencies using the Polygamma functions. 
One can get the result in a finite form, but for our purposes is more useful a low temperature 
expansion (in the ratio T/Aq). In fact, we will be interested in analyzing the line of second 
order transitions starting from T = 0. As anticipated in the Introduction we expect this line 
to become first order at a temperature of about 0.04Ao ^1], therefore an expansion up to the 
order (T/Ao) 2 will be sufficient. 

Notice that for any crystalline structure, the coefficient a depends only on a single q vector. 
Therefore it has a universal structure. In particular the second order points are determined 
by the equation a = 0, and the optimal choice of the vector q along a second order line is 
determined by the condition 

£ =«• (16) 

Cq minimum 

This condition determines |q| leaving the direction of the vector arbitrary. In particular, at 
T = 0, one gets 

qv F ^l.25n2, Sn 2 ~ 0.754A , (17) 

where 6fj,% is the second order transition point and Ao is the BCS gap. Prom this it follows that 
all the vectors appearing in the quantities J and K are of the same length, therefore due to the 
condition of momentum conservation, only a few configurations may appear. In the actual case 
of two antipodal plane waves and for the J integral one has (with q a = — q;, and |q n | = |q&| = q) 

m 

Jo = </(qa,q a ,q a ,qa), Ji = ^(q a ,q a ,qb,qb)- (18) 
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It results 



[3 = (2J + 4J!). 

gp 



(19) 



In the case of K we have three structures 



K = if(q a , q a , q a , q a , q a , q a ), Ki = K(q a , q a , q a , q a , q b , q b ), 

K 2 = if (q a , q a , q b , qb, qb, qb), (20) 



with 



7 = (21^0 + 12^1 + 6^2). 

9P 



(21) 



We are now in the position of evaluating all the quantities of interest at the order (T/Ao) 2 . 
About this point, one can notice that all our expressions are even in T and therefore we are 
really neglecting terms of order (T/Aq) 4 . 



Introducing the re-scaled variables: 



Q 



8/j, 



T -► -r-, Q = qv F -> — , Sfi . 



we get the following results 



5(i 



a(6fi,Q,T) = -l + ^\n 



Q - S/j, 



im 1 



7T 



2 4(Q 2 - V) 3(V - Q 2 ) ' 



j ,*„ OT] _ 9P 1 (, , vr 2 (3^ 2 + Q 2 ) 2 \ 
J (^, Q,T) - - j^-^ 2 [I + _ g2)2 T j , 



(22) 

(23) 
(24) 



Ji(8fi,Q,T) 



99 

K Q {8n,Q,T) = 



1 

Q 

99 



In 



<V + Q 



89- Q 
1 



vr 2 / 2(2^ 2 -Q 2 ) 1 

3^ v (tp 2 - Q 2 ) 2 $9Q n 



Sfi + Q 



89- Q 



64 (8n 2 - Q 2 ) 3 



2 , _ 2 , 27r 2 (5<^ 4 + 10fy 2 Q 2 + Q 4 ) fT1 2 



3V + Q 2 + 



(<J/i 2 - Q 2 ) 2 



(25) 
(26) 



K 1 (8 Ll ,Q,T)=K 2 (8(i,Q,T) 



'2vr 2 (12^ 6 - 14<5// 4 Q 2 + 11<^ 2 Q 4 - 3Q 6 ) vr 2 
+ I — — o ^ttti : o - In 



f 28li 2 - Q 2 1 <5^ + Q 
64<5/x 2 |_(<5// 2 - Q 2 ) 2 28fiQ n 8fi-Q 

8fj, + Q 



Z8n 2 (8n 2 - Q 2 ) 4 



8(i — Q 



(27) 



We have now the expansion of the grand potential up to terms of order A 6 . However for 
the following discussion one would need also the next term in the expansion. Since this is not 
available we will assume it to be strictly positive and we will discuss the results as functions of 
the coefficient 8 of this term. Therefore, our grand potential will be 



- = aA 2 + ^A 4 + ^A 6 + -A 8 , 



Q 

p 4 6 

with a, (3 and 7 from the calculations of this Section and 8 assumed strictly positive. 



(28) 
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3 The general phase diagram 



To study the structure of the phase space let us define a dimensionless grand potential: 

~ n / 6 c 2 1 ,\ 

n= imffi =x \ a+ 4 x+ ** + *V' (29) 



where 



Aq is the BCS gap and 



a= 5Af' 6= 5Af' C= Mf' (30) 
A ^ 2 



The minima of the grand potential are given by the equation 



. (2a + bx + cx 2 + x d ). (32) 
oA Ao 

Notice that the structure of the minima is not changed by dividing by the parameter 5 since it 
is assumed to be strictly positive. Therefore our phase space is a three-dimensional one. We 
want to determine the regions of this space corresponding to first and second order transitions. 
Their intersection will fix the tricritical lines. We start analyzing the roots of eq. (|32|) different 
from zero. The type of solutions with their sign is given in Table ^ Keeping in mind that we 
have the constraint x > for not zero roots, we can easily determine the kind of symmetry in 
each octant, and the result is shown in Fig. ^ 

In particular we see that for a < the symmetry is always broken. This is confirmed also 
by looking at the second derivative of O: 

Aqt^2 = 2 a + 3bx + 5cx 2 + 7x 3 . (33) 

The point x = for a < is always a maximum, whereas for a > is a minimum, and we have 
to decide which is the absolute minimum. If we evaluate the second derivative at a root x r ^ 
we find 



-12a- 4bx r -2cx 2 r , (34) 



which is negative in the octant 1 (a, 6, c > 0). Therefore in this octant the true minimum is at 
x = and we are in the symmetric phase as seen from the analysis in Tabled To decide what 
happens in the octants 4, 5 and 8, the ones with mixed symmetry, it is enough to look at the first 
order points. These are the points where, for a / 0, we have a change in the symmetry. Since 
a ^ 0, the only possible minimum in zero is the trivial one, and these points are determined by 
a change in the number of the real solutions of the cubic 

Cl b C n 1 n , % 

- = a + -x+ -x l + -x A = 0, (35) 
x 4 6 8 

as already evident from Tabled Therefore the first order surface is made up by the points of the 
regions 4,5 and 8 where the discriminant of the cubic equation Q35|) is zero. Before discussing this 
point let us notice that the second order transitions are characterized by a = 0, but not all this 
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a 


b 


c 


Solutions 


Phases 


1 


+ 


+ 


+ 


3 negative 
1 negative, 2 complex 


symmetric 
symmetric 


2 




+ 


+ 


2 negative, 1 positive 
2 complex, 1 positive 


broken 
broken 


3 




+ 




3 positive 
2 complex, 1 positive 


broken 
broken 


4 


+ 


+ 




2 positive, 1 negative 
2 complex, 1 negative 


broken 
symmetric 


5 


+ 




+ 


2 positive, 1 negative 
2 complex, 1 negative 


broken 
symmetric 


6 






+ 


2 negative, 1 positive 
2 complex, 1 positive 


broken 
broken 


7 








2 negative, 1 positive 
2 complex, 1 positive 


broken 
broken 


8 


+ 






2 positive, 1 negative 
2 complex, 1 negative 


broken 
symmetric 



Table 1: The table shows the character of the solutions of the cubic according to the sign of the 
coefficients. Also shown are the phases corresponding to the various type of solutions. 

plane is a second order surface. Regions that correspond certainly to second order transitions 
is the part of the plane o = separating region 1 from region 2. In the other cases, that is 3-4, 
5-6 and 7-8 one has first to discuss the region of the first order transitions. 
To discuss the solutions of (|35|) we bring it to the normal form 

y 3 +py + g = 0, (36) 

where 

y = J + * (37) 

and 

16 o 128 o 8 

p = 2b-— c 2 , q = —c 3 - -be + 8a. (38) 
F 27 ' H 729 9 v ; 

The discriminant of the cubic is then given by 

D = — + — (39) 
4 27 v ; 

and we have one real and two complex solutions for D > and three real solutions for D < 0. 
Therefore for D < we are typically in the broken phase. The situation is made clear in Figure 
121 where we show the curves D = in the regions 4, 5 and 7 for different values of b. 

We can now identify the tricritical lines noticing that they are necessarily at the intersection 
of the plane a = with the surface D = passing through the regions 4, 5 and 7. Evaluating 
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D at a = we get: 

= 5*)- (40) 

The line a = and 6 = for c > is certainly a tricritical line, since it belongs to both to the 
first order and the second order surfaces. As far as the line 6 = 2c 2 / 9 it can be critical only for 
6 > and c < (region 4). In fact, in this case the first order lines touch the surface a = 0, 
as shown in Fig. [21 whereas in the regions 5 and 7 this does not happen (see again Fig. [2J. 
Summarizing the tricritical lines are defined by 

a = p = 0, 7>0 (41) 

and 

96 - 2c 2 = 0, 6 > 0, c < 0. (42) 

The first and second order critical surfaces with their intersections are shown in the three 
dimensional plot of Fig. E3 



4 The second tricritical point for the two antipodal plane waves 

We have now all the elements to analyze what happens in the case of two antipodal plane waves 
along the second order transition line starting from T = 0. These points, as we have discussed 
in the previous Section are defined by the condition 

a(Sfi, T, Q) = 0, (43) 

supplemented by the condition of optimality for Q which amounts to require the stationarity of 
the grand potential f2 

dQ K ' 

We can solve these two equations finding the values of Q along the second order transition line 
and a relation between T and 5fi. Of course, the result is the same as in the case of a single 
plane wave, since as noticed before, the expression of a does not depend on the number P of 
plane waves. The transition line in the plane (5/j,, T) is given for small values of T in Fig. 21 
What is different now from the case of the single plane wave is the behavior of the coefficients (3 
and 7. In fact, looking at Fig. [31 for the single plane wave one has a path going from C (T = 0) 
to A (the tricritical point) staying along the second order surface, whereas in the actual case 
one moves from the region (3 > and 7 > to the region (3 > and 7 < where the second 
tricritical line (see eq. (US) ) 

9/35 - 2 7 2 = (45) 

is met. This is shown in Fig. [3 where the behavior of (3 and 7 along the second order transition 
line is shown as a function of the temperature. In particular we see that for 0.045Ao < T < 0.071, 
(3 is positive and 7 is negative. Therefore in this region another tricritical point is found. The 
transition becomes first order and stays like that till the first tricritical point A. The location 
of the second tricritical point is given by equation (|45|) and the results are illustrated in Fig. H3 
This Figure shows that the location of the tricritical point does not depend very much on the 
dimensionless parameter <5Aq. In fact when it varies between zero and one, the temperature of 
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the second tricritical point varies of about 25%, whereas the value of Sfj, changes of about 1.6%. 
The optimal value of Q stays practically constant in this range. 

We can compare our results with the one obtained in |14j . In this paper it has been found, 
by using the method of the quasi-classical Green's function, that the second tricritical point for 
the two antipodal plane waves is located at T « 0.075 T c » 0.043 Aq. By looking at Figure H3 
this suggests that the actual value of 5 at small temperature is rather small, but this conclusion 
should be confirmed by an explicit evaluation of this parameter. Also it should be noticed that 
the analysis made here using the Ginzburg-Landau expansion along a second order line puts 
this result in a more firm basis, and clarifies on a general ground the existence of a second order 
transition point. 

5 Conclusions 

In this paper we have discussed the phase structure of the LOFF phase in the configuration of 
two antipodal plane waves at low temperature. The interest in this study is because it is known 
that at T = the system undergoes a second order phase transition at a critical value of the 
chemical potential separation 5fi of the two pairing species of fermions. On the other hand the 
transition is first order at the tricritical point Ptric- Here we clarify the reasons for this behavior 
as due to the existence of a second tricritical point. The most important point of the paper is 
that, in the case of a Ginzburg-Landau expansion up to the eighth order in the gap, the phase 
space is essentially three-dimensional and that two tricritical lines exist. In fact we have studied 
in full generality the critical points in this space, under the assumption that the eighth term 
in the Ginzburg-Landau expansion is strictly positive. The existence of the second tricritical 
point is then a simple consequence of the general structure of the phase space when the eighth 
order term becomes important, and it comes about since starting from T = and increasing the 
temperature along the second order transition line, one follows a path in the physical phase space 
which crosses the second tricritical line. The study made in this paper, although motivated by 
a particular situation (the two antipodal plane wave of the LOFF phase) , is in fact much more 
general and it is conceivable that it can be applied to other interesting physical situations. 
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Figure 1: The octants of the space (a, b, c) are illustrated in the two planes corresponding to 
b > (left panel) and b < (right panel). Also the phases in the various regions are indicated. 
S is the symmetric phase, whereas B is the broken one. 




Figure 2: We show the points where D = in the region 4, a > 0, 6 > 0, c < (left panel) and 
in the regions 5 and 7 , a > 0, b < (right panel). In the part of the octants below these curves 
we have D > 0, therefore there are 2 complex and one negative root (see Table^j, and we are 
in the symmetric phase. By continuity the region above the curves D = belongs to the broken 
phase. 
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Figure 3: This three dimensional plot shows the surface of first order points in the relevant 
regions 4, 5 and 7 and the surface of the second order points which lies in the plane a > for 
b and c positive and in the part of the plan b > and c < delimited by the surface of the 
first order points. Also shown are the tricritical lines. In this figure we show also roughly the 
path along the critical line starting at the tricritical point (A) of the LOFF phase down to zero 
temperature, point C. In the path the other tricritical point B is found. The rest of the path 
between B and C is along the second order surface. The arrows in the boxes containing broken 
and symmetric phases indicate that these two phases are respectively behind and in front of the 
critical surfaces. 
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Figure 4: This figure shows the second order transition line in the plane (J/i/Ao, T/Aq). 
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Figure 6: The figure shows the variation of the position of the second tricritical point obtained 
by solving eq. \4°~\ ) inT (left-upper panel) andSfi (right-upper panel) . Also shown is the optimal 
value for the wave vector Q at the tricritical point vs. 5 (lower panel). 
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